The order of the Roberge- Weiss endpoint (finite size transition) in QCD. 
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We consider the endpoint of the Roberge- Weiss (RW) first order transition line present for imagi- 
nary baryon chemical potentials. We remark that it coincides with the finite size transition relevant 
in the context of large Nc QCD and study its order in the theory with two degenerate flavors. The 
RW endpoint is first order in the limit of large and small quark masses, while it weakens for inter- 
mediate masses where it is likely in the Ising 3d universality class. Phenomenological implications 
and further speculations about the QCD phase diagram are discussed. 
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I. INTRODUCTION 

The determination of the QCD phase diagram at finite 
temperature T and baryon chemical /is is one of the out- 
most open problems within the Standard Model of Par- 
ticle Physics. At /.ts = 0, lattice QCD simulations have 
shown the presence of a finite T deconfinement/chiral 
symmetry restoring transition: in the limit of zero or 
infinite quark masses it is associated to a change in the 
realization of some exact symmetry (chiral or center sym- 
metry respectively), hence it must be a real phase tran- 
sition. For finite quark masses no further exact symme- 
tries are presently known, hence the transition may also 
be a smooth crossover. At /i^ ^ the deconfinement line 
starting from fiB = may merge with an analogous tran- 
sition line starting from the T — axis: the latter could 
be first order and a critical endpoint may be present if 
the fiB = transition is a crossover. Large efforts are 
dedicated to the theoretical and experimental search of 
this possible endpoint: unfortunately lattice QCD simu- 
lations at nonzero /is are hindered by the sign problem 
and a direct numerical investigation is hardly feasible. 

A way to partially overcome the sign problem is to 
consider an imaginary chemical potential, fj,B = "ifJ-i- nu- 
merical simulations are feasible and information about 
real hb can be recovered by analytic continuation tech- 
niques p]-[T^|. Definite answers can be obtained regard- 
ing the structure of the phase diagram in the T-fij plane: 
those answers do not give direct information on the finite 
density phase diagram but, as remarked in recent litera- 
ture [ill and as further stressed in the present study, 
they may be relevant to physics at zero or small real /i b ■ 

One of the main features of the T-fij diagram are the 
first order lines met at high T and fixed periodic values of 
m known as Roberge- Weiss (RW) transitions The 
end of such lines is a phase transition itself when mov- 
ing in T at fixed /i/: its order may be relevant to QCD 
phenomenology and is the subject of our study. 

The fiB — partition function in the lattice formula- 
tion is 



where U stands for gauge link variables, Sg is the pure 
gauge action and M is the fermionic matrix: periodic (an- 
tiperiodic) boundary conditions are understood for boson 
(fermionic) fields in the Euclidean time direction. More 
fermion determinants or powers of them may be needed 
depending on the flavor spectrum. 

In the quenched limit (no fermion determinant) the 
action is symmetric under multiplication of all temporal 
links at a given time by an element of the center of the 
gauge group Zn, = {e'^W^fc ^ ^ = 0,...iVc - l}. This is 
known as center symmetry and gets spontaneously bro- 
ken at the deconfinement transition. The Polyakov loop 
L is not invariant and serves as an exact order parame- 
ter: (L) is nonzero and proportional to a center element 
in the deconfined phase. The fermion determinant breaks 
center symmetry explicitely: {L) is always nonzero and 
real, but it still rapidly increases at the transition. 

In order to discuss the phase diagram in T-fj,j plane, 
we introduce the dimensionless variable 9q = Im(/Xg)/T, 
where /iq = /i^/S is the quark chemical potential: /i/ can 
be interpreted as a constant U{1) background field in the 
Euclidean time direction or, equivalently, as a twist in the 
fermionic temporal boundary conditions by an angle 9q, 



Z{T,0q)= VUe-^'^^'^^detM[U,0q] 



Z{T)= / VUe-^^'^^^deiMp] 



(1) 



(2) 



A transformation 9q ^ 6q + 2Trk/Nc can be exactly can- 
celled by a center transformation leaving both Sg and 
the functional integration invariant, hence the free en- 
ergy is periodic in 9q with period 27r/A^c instead of 
the expected 27r. Such periodicity is smoothly realized 
at low T, while at high T phase transitions occur for 
9q ~ {2k + 1)tt/Nc and k integer, at which (L) jumps 
from one center sector to the other: the phase of L and 
other physical observables (e.g. the baryon density) are 
discontinuous at these points The emerging picture 
for the T-9q phase diagram is that of a periodic repetition 
of first order lines (RW lines) in the high T regime, which 
must disappear in the low T regime, hence they have an 
endpoint at some temperature Trw. There is numeri- 
cal evidence that the analytic continuation to imaginary 
/is of the deconfining/chiral restoring line (present for 
real /is) touches the RW line right on its endpoint and 



2 



then repeats periodically in the T-0q plane: we shall com- 
ment on this issue later. The deconfining temperature in- 
creases with Hi, hence Trw > Tc, where Tc is the critical 
temperature at fiB =0. 

RW lines corresponds to points where a Z2 symmetry is 
spontaneously broken: two center sectors are equivalent 
but one of them is selected. Let us consider in particular 
the 9q = -K line. As for 9q = there is symmetry with re- 
spect to complex conjugation of link variables, i.e. charge 
conjugation, however in the high T phase for 9q = n the 
system selects one of two center sectors in which (L) is 
complex and charge conjugation is spontaneously broken: 
Tptw is the critical temperature at which charge symme- 
try breaks, hence it corresponds to a real phase transition 
independently of the quark spectrum. The RW line is 
analogous to the first order line lying along the T axis of 
an Ising system: the endpoint plays the role of the Ising 
critical T and {9q — tt) that of the magnetic field. It is 
worth stressing that instead at 6*^ = (L) is always real 
and charge conjugation stays unbroken. 

The RW endpoint can be given a different interpre- 
tation. Setting = TT is like switching from antiperi- 
odic to periodic boundary conditions for fermion fields in 
the Euclidean time direction. The system can be viewed 
as a usual thermal system in presence of an imaginary 
chemical potential or, naming differently the axes , as a 
T = system with one spatial dimension compactified: 
in the second case the RW endpoint defines a critical 
size Lc = 1 /Trw of the compact dimension, below which 
charge conjugation is spontaneously broken. 

Such transition has been studied in recent litera- 
ture [l8l - [2l| and is relevant to large Nc QCD, in par- 
ticular to large Nc orientifold planar equivalence [22] . 
That states that QCD with fermions in the antisym- 
metric representation, QCD(AS), which coincides with 
ordinary QCD for Nc = 3, is equivalent, in the large Nc 
limit and in the charge-even sector of the theory, to QCD 
with fermions in the adjoint representation. Such equiva- 
lence is guaranteed if the charge conjugation symmetry is 
not spontaneously broken in QCD(AS) (actuall y th e 
condition may be not so strict, as discussed in Ref. j24|). 
That fails in presence of a compact dimension, below a 
critical compactification radius, which is nothing but a 
different mapping of the endpoint temperature Trw 

The existence of a phase transition at the RW end- 
point is fixed by symmetry, but its nature is not. If it is 
second order, symmetry suggests that it belongs to the 
universality class of the 3d Ising model; the correspond- 
ing critical behaviour could in principle have influence 
also far from the endpoint, e.g. for zero density QCD 
right above Tc, as suggested by recent literature [lOlllSj. 
However it may also be first order: this is the case in the 
infinite quark mass limit, where the transition coincides 
with the pure gauge SU(3) thermal transition. In general 
the answer may depend on the flavor spectrum. 

In the first order case, the RW endpoint is actually 
a triple point, with two further first order lines depart- 
ing from it for non-zero values of the magnetic variable 



(think of the 3d 3-state Potts model as a similar exam- 
ple). It would be natural to identify this departing line 
with (part of) the analytic continuation to imaginary fj,B 
of the physical line, thus explaining why it meets the RW 
line right on its endpoint. The departing line could reach 
the fij — axis or have a second order critical endpoint 
arbitrarily close to it, thus with a strong influence on the 
properties of strong interactions right above Tc. 

In the following we show results regarding QCD with 
two degenerate flavors. Earlier lattice results [S] on small 
lattices suggested a second order nature, supp orted re- 
cently by the analysis of effective models [15]. Instead 
we provide evidence that the transition is first order for 
low enough masses and that it weakens for intermediate 
quark masses, where it could be second order. 

II. NUMERICAL RESULTS 

We have investigated QCD with two degenerate fla- 
vors, adopting the standard plaquette action, the stan- 
dard staggered fermion formulation, and using a Rational 
Hybrid Monte Carlo algorithm. Two values of the quark 
mass have been explored, aniq = 0.075 and ama = 0.025, 
the latter coinciding with that used in Ref. Q. In or- 
der to perform a finite size scaling analysis at the critical 
endpoint, we have made simulations on lattices x Lt 
with Lt =4 and = 8, 12, 16, 20, 32. We have worked 
at fixed 9q = t: and the temperature T = 1/ {Lta{p, mq)) 
has been changed by tuning the inverse gauge coupling /?. 
Collected statistics are of the order of 50 — lOOK trajecto- 
ries for the /? values closest to the critical point. Simula- 
tions have been performed on a computer farm in Genoa, 
apart from those on the biggest lattice, Lg ~ 32, for 
which the apeNEXT facilities in Rome have been used. 
Both resources have been provided by INFN. 

Since we work at 9q = tt, we have chosen the imaginary 
part of the Polyakov loop as an order parameter: it is 
not invariant under charge conjugation and develops a 
nonzero expectation value for T > Tr^v (another possible 
choice could have been the imaginary part of the baryon 
density [lEIl^): actually we consider its modulus for the 
finite size scaling analysis, similarly to what is done for 
the magnetization of an Ising system. Its susceptibility 

X^£H(Mi)^)-(|Im(L)|)2) (3) 

where L is the spatially averaged Polyakov loop trace 
(normalized to Nc), is expected to scale as follows: 

X = L]^" HrLVn ■ (4) 

That means that the quantities x/^s^"^ measured on dif- 
ferent lattice sizes should fall on the same curve when 
plotted against tL^'^ . The critical indexes are v ^ 0.63 
and 7 ~ 1.24 for the 3d Ising model, while a first or- 
der transition is effectively described, in three spatial 
dimension, hy v = 1/3 and 7 = 1. In the following 
(/? — /3rw) will be used in place of t = (T — Tp{vv)/7rw: 
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FIG. 1: Scaling of the reweighted susceptibility of the imaginary part of the Polyakov loop for am = 0.025 (left column) and 
am = 0.075 (right column) according to the Ising 3d universality class (upper row) or to first order critical indexes (lower row). 




FIG. 2: Reweighted plaquette distribution at /3c as a function of Lg for am — 0.025 (left) and am = 0.075 (right). 



that can be done close enough to the critical point. In 
Fig. [T] we show how these two scahng ansatze work for 
the two masses explored; the plotted susceptibilities have 
been obtained by Ferrenberg-Swendsen reweighting. At 
the lowest quark mass the scaling with 3d Ising critical 



indexes clearly fails (upper-left figure), while first order 
scaling is good (lower-left figure), in both cases we have 
used /3rw = 5.33885: these results are in disagreement 
with those of Ref. \3\ and this is likely due to the small 
volumes Lg = 6, 8 and to the non exact molecular dy- 
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namics algorithm used in Ref. Q. 

Instead at the larger quark mass, arrig = 0.075, first 
order scaling is not good (lower-right figure) while better, 
even if not perfect, agreement is found with the 3d Ising 
critical behaviour (upper-right figure). In this case the 
critical coupling has been set to /3rw = 5.3965. 

Consistent results are obtained for the plaquette dis- 
tribution at the critical couplings, reported in Fig. [21 For 
airiq = 0.025 a double peak structure develops, becoming 
sharper and sharper as the volume increases; a similar be- 
haviour (not shown in the figure) is observed for the chiral 
condensate, which is also discontinuos at the transition. 
This is not the case for am, = 0.075: we cannot exclude a 
very weak first order transition with discontinuities bet- 
ter visible on volumes larger than those explored by us, 
but we conclude that at amq ~ 0.075 the transition is 
surely much weaker than at amq = 0.025. 

The simplest scenario compatible with our results is 
that the RW endpoint is first order both for very heavy 
or light quarks, while it weakens and could be second 
order for intermediate masses. Note that aniq = 0.025 
corresponds to a pion mass already quite larger than the 
physical one, therefore we predict a first order RW end- 
point also for physical quark masses. Our results should 
be checked closer to the continuum limit using improved 
actions and/or larger values of Lt; also a check with dif- 
ferent fermion formulations would be important. Nev- 
ertheless, if this scenario is confirmed, it has important 
phenomenological consequences and gives rise to further 
speculations that we discuss in the next Section. 



III. DISCUSSION AND SPECULATIONS 

A first order RW endpoint implies further first order 
lines departing from it. One of those lines is part of the 
analytic continuation of the physical transition line: the 
fate of this first order line is of great importance. Two 
possibilities can be realized: 

i) It ends in a critical endpoint before reaching the 
Hi = axis, as sketched in Fig. |3l However the endpoint 
could be close to the axis and thus have strong influence 
on zero density physics slightly above T^, 

ii) It reaches the /z/ = axis and possibly goes through 
it, implying a first order transition at zero density and 
also for small real chemical potentials. 

It is reasonable to assume that the extension of the first 
order line increases as the strength of the RW endpoint 
transition increases. Hence, according to our results, the 
second possibility should be more likely for very heavy or 
very light quarks: indeed the possible presence of a first 
order transition in the chiral limit of Nf = 2 QCD is not 
excluded by recent lattice studies [1^, [2^ . 

The strengthening of the RW endpoint transition for 
large quark masses is understandable, since the symme- 
try which is spontaneously broken at the endpoint is 
a remnant Z2 subgroup of the original center symme- 
try, which is exact in the quenched limit. Instead the 




FIG. 3: A possible sketch of the T — jj.^ diagram in which 
the first order line departing from the RW endpoint ends at 
a critical temperature before reaching the fj,^ = axis. 



strengthening in the chiral limit is an unexpected phe- 
nomenon, which is likely linked to the interplay with chi- 
ral degrees of freedom and which should be investigated 
by future studies. It is reasonable to assume that such 
phenomenon takes place also for different number of fla- 
vors: then further speculations follow. 
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FIG. 4: Sketch of a speculative phase diagram for Nf = 3 
case, inspired by our results for Nf — 2. Thick solid lines 
represent the border of first order transition surfaces. 

Consider the case of three degenerate flavors, Nf — 3, 
where it is known that, at /^s = 0, the transition is first 
order for small or large quark masses, with two critical 
quark masses rric and delimiting an intermediate re- 
gion where a crossover is present (think of the Columbia 
plot). It is possible to interpret the two first order re- 
gions in terms of the realization of possibility ii) described 
above: for m < rric or m > m'^ the RW endpoint is 
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so strong that the departing first order line reaches the 
Hb = axis; the same does not happen instead for in- 
termediate quark masses rric < m < m'^. This conjecture 
is sketched in Fig. 21 representing a speculative phase 
structure for Nf — 3 in the T-m-fj,^ diagram (imaginary 
chemical potentials corresponds to /i^ < 0): a first or- 
der surface departs from the line of RW endpoints and 
extends enough to reach = only for low or high m. 

The diagram in Fig. |4] is in contrast with a standard 
scenario in which the first order region seen for /i^ = 
and m < nic is the intersection with a first order surface 
related to high /i^, low T physics and to which the criti- 
cal endpoint of the QCD phase diagram belongs, so that 
the critical endpoint is on the same line as the m — rric, 
/x^ = point. In Fig. |4] the high /i^, low T surface is plot- 
ted as well (this is a simplification, since more transition 
surfaces may be present at high /i^), but it is not related 
to the first order transition seen for /i^ = 0, to < TOc, so 
that if a QCD critical endpoint exists, it is not related to 
the chiral critical surface. 

Fig. 2] is a conjecture for Nf = 3 inspired by our 
results for the Nf = 2 RW endpoint: is it a reasonable 
guess? Future numerical studies can clarify that, at 
least for the /i^ < side. Moreover, if we consider how 
TOc changes as fi^ is slightly increased from zero, then 
according to Fig. 2] it should decrease, while according to 
the standard scenario it should increase: the authors of 
Rcfs. [27] have provided evidence showing that TOc(//^) 
is indeed a decreasing function of /i^, at least for small 
/i^, thus supporting a scenario like that shown in Fig. 21 

To summarize, we have provided evidence that the end- 
point of the RW transition for Nf = 2 is first order for 



small or high quark masses, while it weakens and could 
be second order for intermediate masses. A first order 
RW endpoint implies two further first order lines depart- 
ing from it: it is natural to identify those lines with part 
of the analytic continuation of the physical line, thus ex- 
plaining why the latter and the RW line are connected to 
each other. The range of light quark masses, for which 
the RW endpoint is first order, includes physical quark 
masses, it is therefore of great importance to understand 
what is the fate of the further first order line departing 
from the endpoint: it may reach the zero density axis 
or have a critical endpoint arbitrarily close to it, which 
could have great infiuence on the physics of strongly in- 
teracting matter right above the deconfinement transi- 
tion. We will clarify this point in the future, by studying 
how fast the first order observed at the RW endpoint 
weakens as /i/ moves towards the /i = axis. A careful 
check of our results closer to the continuum limit should 
be performed as well. Assuming that similar results may 
be found for different flavor numbers, we have made a 
conjecture about the QCD phase diagram which, even 
if rather speculative, can be carefully checked by future 
studies at imaginary fj,B- In particular one should clarify 
if a phase structure like that sketched in Fig. |4l is really 
valid, at least for < 0. 
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